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^ ■ ABSTRACT 

Globular cluster systems (GCS) frequently show a bi-modal distribution of the cluster 
^ . integrated colours. This work explores the arguments to support the idea that the same 

' feature is shared by the diffuse stellar population of the galaxy they are associated 

I with. The particular case of NGC 1399, one of the dominant central galaxies in the 

■ Fornax cluster, for which a new B surface brightness profile and {B — Rkc) colours are 

presented, is discussed taking advantage of a recently published wide field study of its 
GCS. The results show that the galaxy brightness profile and colour gradient as well 

■ as the behaviour of the cumulative globular cluster specific frequency, are compatible 

with the presence of two dominant stellar populations, associated with the so called 
"blue" and "red" globular cluster families. These globular families are characterized 

pi i' by different intrinsic specific frequencies (defined in terms of each stellar population) : 

Sn = 3.3 ± 0.3 in the case of the red globulars and Sn = 14.3 ± 2.5 for the blue ones. 
We stress that this result is not necessarily confiicting with recent works that point 
out a clear difference between the metallicity distribution of (resolved) halo stars and 
^ , globulars when comparing their number statistics. The region within 0.5 arcmin of 

the centre shows a deviation from the model profile (both in surface brightness and 
colour) that may be explained in terms of the presence of a bulge-like high metallicity 
, component. Otherwise, the model gives an excellent fit up to 12 arcmin (or 66.5 Kpc) 

^ • from the centre, the galactocentric limit of our blue brightness profile. The inferred 

\ specific frequencies imply that, in terms of their associated stellar populations, the 

formation of the blue globulars took place with an efficiency about 6 times higher than 
that corresponding to their red counterparts. The similarity of the spatial distribution 
of the blue globulars with that inferred for dark matter, as well as with that of the 
X ray emiting hot gas associated with NGC 1399, is emphasized. The impact of a 
relatively unconspicuous low metallicity population, that shares the properties of the 
blue globulars, as a possible source of chemical enrichment early in the formation 
history of the galaxy is also brieffy discussed. 
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1 INTRODUCTION 

A large volume of arguments supports the idea that the 
current structure of most galaxies has been determined by 
events that occurred at the very early epo chs of the Uni- 
verse (see, for example, iBaueh et alJlTooM and references 
therein) . In this context, and being the oldest stellar sys- 
tems for which a reliable chronology exists (at least in the 
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Milky Way) , globular clusters (GC) have become the subject 
of an increasing amount of research that aims at identify- 
ing the link between their formation and that of their host 
galaxies themselves. 

Even though the fact that GC may be used as probes 
of th e formation of a ga laxy has been known for a long time 
(e.g. lEggen et aLlligG^l and that the knowledge of globular 
cluster systems (GCS), in our and in other galaxies, has been 
steadily growing, a unified picture that explains both the 
small scale phenomena that lead to the cluster formation, 
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and also their overall properties as systems, is still missing. 

Recent reviews give good summaries of the situa tion 
fe.g. lHarri£l200.?l : lKissler-Patiel20 02:'Kissler-Patig'2003"l and 
compare the features of several astrophysical scenarios that 
invoke "external" or "internal" phenomena that provide a 
reasonable confrontation with observations. It seems clear, 
however, that some aspects of the proble m require a cos- 
molog ical frame as already suggested by IPeebles fc Dickj 
ilOGS") . More recently, other atte mpts assume hierarchical 
CDM models iBeaslev et al .l2002h that appear as a promiss- 
ing approach to further incorporate the effects of different 
and co-existing (rather than excluding) processes that end 
up with the formation of a GCS as, for example, the role of 
GC as a possible source of the re-ionization of the Universe 
jRicottil200i ISantojl2003l ) . 

All those scenarios, however, must take into account 
some relatively well established facts as boundary condi- 
tions, namely: 

- Stellar clusters that fulfill the properties expected for 
a "young" globular cluster are currently formed in violent 
events, such as g alaxy mergers, as originally pointed out by 
ISchweizeil (Il987l) . These kind of phenomena, on a variety of 
scales, were probably more frequent during the early phases 
of galaxy formation. 

- GC formation seems connected with major stellar for- 
mation events. This fact is well supported by the rela- 
tively narrow range of the ratio of the number of GC to 
galaxy luminosity, as well as by the co rrelation of cluster 
formation with high star formation rate iLarsen fc Richtleil 
I2OOCII) . This ratio can be quantified, for example, using the 
so called "specific globular clu ster frequency", Sn, following 
iHarris fc van der Berghl (^^). We point out that a revision 
of several archetypical high Sn galcixies have lowered previ- 
ous estimates of this parameter by half, or eve n less. (e.g. 
lOstrov. Forte fc Geisleilll998l: |^rtee£al]|2Qo2) . 

- The Milky Way GC metallicity distribution shows a 
bi- modal structure. This feature seems shared by other 
galaxies that show bi- modal colour distributions (e.g. 
iGeisler. Lee fc Kimlll996l) . Although the well known colour- 
age- metallicity degeneracy prevents a stronger statement 
about the significance of the colour distribution, the assump- 
tion of very old ages for GC points to metallicity as the main 
driver of the colour distribution. 

- In the case of NGC 1399, the so called "red" GC exhibit 
a distinct behavior, when compared with the "blue" GC, 
both in terms of their spatial dis tribution and kinematic 
behaviour fe.g. lRichtler et alJliool : R2004 in what follows). 

In this picture, the connection between the GCS and the 
stellar pop ulation in a galaxy is not com pletely clear. Some 
differences fForte. Strom fc Strom| [l98l[) and some similar- 
ities (Forbes fc Forte 2001) have been pointed out over the 
years. Th is last paper, in pa rticular, reinforces previous ar- 
guments iForbes et al.lll997r) that suggest a connection be- 
tween the red GC and the galaxy halo stars based on the 
similarity of the red GC colours with those observed in the 
inner regions of a galaxy. A si gnificant contribution to the 
subject has been presented bv lHarris fc Harris i2002h (and 
references therein) who discuss the characteristics of the re- 
solved stellar population of NGC 5128 and present a com- 
parison with the GC in that galaxy. These authors also dis- 



cuss a possible connection between GC specific frequency 
and metallicity. 

This paper explores the arguments to support (or re- 
ject) the idea that, in the bi-modal GCS, the "diffuse" galaxy 
stellar population also exhibits a dual nature reflecting the 
existence of two major stellar populations that share the 
properties of each GC family and are typified by a given 
"intrinsic" Sn, integrated colour, and spatial distribution. 

We emphasize that this work presents an approach that, 
in principle, makes no assumptions about the nature of inte- 
grated colours (in terms of ages or metallicities) and finally 
aims at matching the shape of the observed galaxy surface 
brightness, colour gradient and cumulative Sn (with galac- 
tocentric radius) by properly weighting each of the diffuse 
stellar components. These stellar populations, defined just 
in colour terms, then might include a range of age and/or 
metallicity as well as of GC specific frequencies. 

NGC 1399, the central galaxy in the Fornax cluster, 
appears as a good candidate on which to perform this kind 
of analysis since, as in other luminous ellipticals, its GCS 
is well populated and can provide statistically significant 
results. In addition, it is relatively nearby. 

A number of papers have dealt with the clearly bi- 
modal GC colour distribution in this galaxy taking advan- 
tage of the high metallicity sensitivity of the colour in- 
dex defined in terms of the C and Ti band s of the Wash- 
ington photometric svstem i Canternalll976l) . These early 
works iGeisler fc Fortdll99d : lOstrov. Geisler fc Fort j|l993l: 
iQstrov et al.lll9 98') have a relatively small areal coverage. 
In turn. Tpireli et al. (2003) (D2003 hereafter) increased the 
surveyed area to about 900 sq. arcmin by means of mosaic 
CCD observations and determined GC areal density profiles 
up to a galactocentric radius of about 20 arcmin. These den- 
sity profiles, and the surface photometry (Johnson's B and 
Kron-Cousins R filters) presented in this work are the basis 
of the subsequent analysis. An attempt to disentangle the 
stellar populations in NGC 1399 was presented in D2003 
although, as described later, with a somewhat different ap- 
proach. 

The structure of the paper is a follows: Section |21 de- 
scribes the result of composing two stellar populations and 
the variables that govern the shape of the galaxy surface 
brightness profile as well as the derivation of the intrinsic 
specific frequencies of each GC population; The character- 
istics of the GCS system associated with NGC 1399 and a 
justification of the adopted distance modulus are presented 
in Section |3 a new blue profile and {B — Rkc) colors for 
the galaxy are presented in Sectional The areal density run 
with galactocentric radius for both the red and blue GC, as 
well as the derivation of their volumetric density profiles, are 
discussed in Sections |S] and |H] respectively; Section |7] gives a 
confrontation of the observations with the properties of the 
model in terms of the shape of the surface brightness pro- 
file, colour gradient and cumulative GC specific frequencies 
as a function of galactocentric radius. A caveat about the 
approach is commented on in Section |H| A first approach to 
the expected metallicity distribution function of the galaxy 
halo stars is described in Section |5| The total stellar mass 
as well as the GC formation efficiencies are discussed in Sec- 
tion ll(JI A comparison of the GC spatial distribution with 
that inferred for dark matter (R2004) and hot gas obtained 
from X ray observations are given in Section [TTI A discus- 
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sion of the results and possible implications on the galaxy 
formation scenario are included in Sections 1121 and 1131 re- 
spectively. 



2 A COMPOSITE SURFACE BRIGHTNESS 
PROFILE 

In what follows we assume that both the red and blue glob- 
ulars in NGC 1399 have an associated diffuse stellar popula- 
tion characterized by their own intrinsic specific frequencies. 

A generalisation of the Sn parameter to any photomet- 
ric band, and in particular, for the B (Johnson's blue) band, 
leads to a composite luminosity (in M{B) = — 15 units): 



L{B) = + N{BGC) 



Sb{RGC) Sb{BGC) 



(1) 



The blue band is chosen, in this case, since we will at- 
tempt to match the shape of the B surface brightness pro- 
file presented in Section 4. If N(RGC) and N{BGG) are the 
projected areal densities of the red and blue GC respectively 
and {V — Mv)o is the (interstellar extinction) corrected dis- 
tance modulus, for any galactocentric radius: 



{V - Mv)o = ^13- A{B) - M{B) 



(2) 



where fis is the galaxy surface brigntness, A{B) is the in- 
terstellar extinction along the line of sight in the blue band, 
and M{B) = ~2.5log{L{B) x 10'') is the areal absolute blue 
magnitude, if L{B) is given in absolute brightness units of 
M{B) = —15.0; and if both intrinsic specific frequencies are 
assumed to remain constant with galactocentric radius, then 
it follows: 

^iB = {V- Mv)o + A(B) + 2.Uog[SB{RGG)\ 

-2.5%[iV(i?GC) + ^^^^^] (3) 

where, C{B) is the ratio of the specific frequencies of the 
blue to red GC populations, Sb{BGC)/ Sb{RGG). 

The expression for /is shows that, if the areal densities 
of both cluster families as a function of galactocentric radius 
are determined from observations, the shape of the surface 
brightness profile will depend only on the C{B) parameter. 

Aiming at obtaining the best fitting to the observed 
^B in terms of the globular specific frequencies we then 
proceeded in two steps. First, we looked for the C{B) pa- 
rameter that makes the quantity ^b + 2.5log[N {RGG) + 
N{BGC)/C{B)\ constant with galactocentric radius and, 
from this constant, after adopting a given distance modu- 
lus and interstellar reddenig, we obtain Sb{RGG). In turn, 
Sb{BGG) comes from C{B) and from this last parameter. 

This procedure differs from that given in D2003 in the 
sense that:a) Only the profile shape is fitted (colour gradi- 
ents are predicted and then compared with observations; b) 
Constant intrinsic specific GC frequencies are assumed (i.e., 
not forcing a particular functional depedence with galac- 
tocentric radius) ;c) The derived surface brightness is con- 
fronted with an observed profile that relies on genuine sky 
brightness measurements (as opposed to the adoption of a 
preferred value as in D2003). 

In the following sections we describe all the terms in- 
volved in the ^b expression that will be used to fit the ob- 
served blue profile. 




Figure 1. {C — T\) colour distribution for 1246 globular cluster 
candidates associated with NGC 1399 and within a galactocentric 
radii of 2.0 to 7.0 arcmin. The sample includes globular candidates 
with Ti (or Rkc) magnitudes between 21.0 and 23.2; colour bins 
are 0.1 mag wide. Two well defined modal colours are seen at 
(C - Ti) = 1.28 and 1.77. 



3 THE NGC 1399 GCS AND ITS DISTANCE 
MODULUS 

3.1 Bi-modality 

The possible existence of a bi-modal (C — T\) colour 
dist ribution in the NGC 1399 GCS was suggested 
bv [ Ost rov, Grgisler fc Fortd (^g^) and confirmed in 
lOstrov et al.l (1199811 bv means of (C - Ti) vs. (Af - Ti) 
two colour diagrams. This last work noted that, except in 
the very inner regions of the galaxy, the modal colours of 
both cluster families remained at constant values indepen- 
dent of galactocentric radius. This feature was also noticed 
by D2003 who extended the survey to about 20 arcmin in ra - 
dius. As noted in this last work, and based on lCeisleil il996f) . 
the Ti and Ruc bands are practically identical and in what 
follows we adopt (C — Ti) as equivalent to (C — Rkc)- 

Figure shows the (C — Ti) colour distribution for a 
sample of 1246 GC candidates contained within an annular 
region between 2 and 7 arcmin of the galaxy centre, taken 
from the Forte et al. (2002) and D2003 data. The annular 
boundaries of this region were chosen in order to maximize 
the number of GC while keeping a low contribution from 
field interlopers. In turn, the magnitude range covers from 
Ti = 21.0, to the position of the GC luminosity function 
turn-over at T\ = 23.2. We ad opt this val ue as a compro- 
mise between those derived bv lCeisler fc Forte (1990) and 
lOstrov et aP (|l998). 

GC candidates brighter than T\ = 21, that were left 
out of the analysis, constitute a small fraction of the whole 
cluster population (less than 9%) and their nature remains 
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to be clarified as th ey exhibit a uni-mo dal colour distribu- 
tion (noted both in lOstrov et aljll998l and D2003) and a 
higher radial velocity dispersion than that characterizing 
their fainter counterparts (from R2004 data). 

A convolution of the original data, that in Figure ap- 
pear binned in 0.10 mag intervals, with a gaussian kernel 
with a dispersion of 0.05 mag (similar to the photometric 
errors) leads to (C — Ti) modal colours of 1.28 and 1.77 for 
the blue and red globulars, respectively. In what follows, and 
consistent with D2003, we adopt a (C — Ti) colour of 1.55 as 
a formal boundary between the blue and red GC (their as- 
sociated diffuse stellar populations will also be referenced as 
"b lue" and "red" in wh at follows). The colour excess maps 
bv lSchlegel et al.l Jl998ll give E{B — V) = 0.015 that trans- 
lates to E{C - Ti) = 0.030 along the line of sight to NGC 
1399. The resulting intrinsic colours are {C — Ti)o = 1.25 for 
the blue GC and (C - Ti)o = 1.74 for the red GC. In order 
to derive a. {B — R) colour gradient and the specific frequen- 
cies of both cluster families in the V band, S„, by means of 
the best profile fitting to the observed hb (Section [SJ, we 
used the following transformations: 

(B-\/)o = 0.508(C-Ti)o + 0.064 (4) 
(B-i?)o = 0.704(C-ri)o + 0.269 (5) 

These last relations were obtained by means of bi- 
sector fits to the intrinsic colours of MW globulars 
iReed. Hesser fc Shawll^l988^ that also give a good repre- 
sentation ofJrilrmsk_£ol(^^ galaxies with redder globu- 
lars (see lForbes fc FortdlioOll) and yield {B - V)o = 0.68, 
(B - R)o = 1.13 for the blue GC and {B ~ V)o = 0.95, 
{B — R)o ~ 1.49 for the red GG with an overall uncertainty 
of ^ ±0.03 mag. 

3.2 Distance modulus 

We derive a distance modulus to NGC 1399 by means of 
the turn-ove r magnitude of t h e blue GC, Bto = 24.32 mag, 
obtained by iGrillmair et al.l 1^9^ on the basis of their 
HST photometr y. In order to avoid cor rections for metal- 
licity effects (e.g lAshman. Conti fc Zepll lT995l we compare 
that magnitude with the mean blue absolute magnitude of 
M W globulars wi th (C — Ti)o colours in the same range. 
The lHarrij il996ri compilation includes 35 globulars with a 
mean (C - Ti)o = 1.15 (i.e. slightly bluer than the N1399 
mean blue GC colours) and a mean absolute magnitude 
M{B) = -7.15 ± 0.3. Adopting A{B) = 0.06 mag (con- 
sistent with the E{B — V) excess mentioned above) leads 
to {V - M^)o = 31.4 (corresponding to 19 Mpc). This dis- 
tance modulus agrees with that adopted by R2004 who give 
a summary of other distance estimates available in the lit- 
erature. 



4 BLUE SURFACE BRIGHTNESS PROFILE 
AND {B Rkc) colours FOR NGC 1399 

The determination of galaxy surface brightness profiles over 
large galactocentric radii faces well known problems con- 
nected with the adoption of a proper sky level and fiat 
fielding (particularly severe in wide field instruments) . Small 
field telescopes, on the other side, are more manageable in 



terms of fiat field quality but, in the case of objects with 
large angular sizes, do not allow the measurement of a mean- 
ingful sky brightness. 

In this work we used the 2.15 m telescope at the Com- 
plejo Astronomico El Leoncito (San Juan, Argentina) in 
combination with a focal reducer that gives an effective field 
of 9 arcmin on the CCD detector. Simultaneous and contin- 
uous measurements of the sky brightness wore obtained by 
means of a zenithal airglow sensor (fScheer igS?! that mea- 
sures the intensity of atmospheric molecular emission lines. 
These measurements were correlated with B and Rkc sky 
brigtness estimates made on CCD frames taken 40 arcmin 
to the north of the galaxy centre (and intertwined with the 
galaxy field images). 

The correlation took into account the differences in air- 
mass between the zenith and the line of sight to the galaxy 
as well as atmospheric extinction (derived during the run: 
fcs = 0.29; fcfl = 0.08). The sky brightness to be subtracted 
from each galaxy field was then obtained by interpolating its 
temporal beahviour derived from the airglow vs. sky frames 
correlation. 

Sky flats were obtained inmediately after sunset and 
provided an overall flat field on the order of 0.5%. 

The observation sequence included sets of four images 
(15 min exposures) centreed on the galaxy and off set frames 
taken along a position angle of 45 degrees (N to E) and 
reaching a galactocentric distance of about 12 arcmin. This 
direction is clean of bright galaxies that may contaminate 
the NGC 1399 outer profile. 

Typically, the B sky magnitude ranged from 23.1 to 
23.3 mag per square arcsec during the night. We note that 
the observing run took place in December 1997 and during 
a period of exceptionally low airglow level coincident with a 
minimum of solar activity. 

Images centered on the galaxy were processed with the 
ELLIPSE routine within IRAF^ and delivered brightness 
profiles that reach 240 arcsec from the galaxy centre. The 
blue profile was extended up to 510 arcsec using three 20 
minute exposure images obtained with the 4m Blanco Tele- 
scope at CTIO (described in Forte ct al. 2002). These last 
images are not suitable for an estimate of the sky level and 
where then matched to our own profil es by means of the 
method described bv lCaon et al.l l|l99l). 

The much higher sky brightness in the R^c band pre- 
vented a similar approach, restricting the limiting galacto- 
centric radius to 240 arcsec with a maximum estimated error 
of ±0.04 mag at that radius. 

In the case of the off-centre blue frames, aperture pho- 
tometry was carried out within 25 square windows (5 arcsec 
on a side) located on a region of the images were the fiat field 
was considered better than 0.2% and far from bright stars 
(in order to avoid spurious scattered light). These measures 
provide the surface brightness at mean galactocentric radii 
of 410, 470, (overlapped with the extended profile), 605 and 
730 arcsec. 

The final composite blue profile is shown in Figure |5| 



^ IRAF is distributed by National Optical Astronomy Observa- 
tories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under co-operative agreement with 
the National Science Foundation. 
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Figure 2. Blue surface brightness profile for NGC 1399 as a 
function of elliptical radius (re = Vofe). The bars correspond to 
the error of the mean magnitude. Errors for galactocentric radii 
smaller than 400 arcsec are comparable to the size of the plotted 
symbols (±0.03 mag.) The solid line is a Sersic profile. The long 
dashed and short dashed lines correspond to one a and to three 
a uncertanties of the sky level. 



and listed in Table 1 where [B — Ruc) colour indices are 
also given up to an elliptical galactocentric radius of 228.8 
arcsec. The elliptical galactocentric radius (re = VaS) was 
computed using the ellipticities delivered by ELLIPSE and 
adopting a constant ellipticity e = (1 — h/a) = 0.13 for the 
two outermost off-set frames. Error bars depict the error of 
the mean magnitude derived on the outer CCD frames. 

Figure |5| also includes, as a reference, the fit of a Sersic 
law, ^iB = bo + bi{re/a)", with bo = 3.74, bi = 1.086, a = 
4.82 X 10"" and n = 0.0878. 

Large field blue su rface brightness profiles for NGC 1399 
have been published bv lKilleen fc BicknelJ il988l) (who com- 
bine t heir obse rvations with V values from lSchombertll98(if l 
and bv lCaon"e t al. ( 1994). A comparison between these pro- 
files shows marked differences that increase with galactocen- 
tric radius. Figure|21shows the brightness residuals as a func- 
tion of re^'^^ when comparing the profiles of those authors 
with ours (that, approximately, falls in the middle of both). 
This diagram also includ es a comparison wit h the V bright- 
ness profile presented bv lHilker et al.l lll999ll . by adopting a 
constant difference (B—V) — 0.93, indicating a good overal 
agreement (although some residual differences might be ex- 
pected if a colour gradient were present). 

D2003 present Rkc and C profiles that reach some 20 
arcmin in galactocentric radius. These profiles were obtained 
by adopting a preferred sky brightness on the basis of mea- 
surements at the edge of the images. In particular, their 7?^^ 
profile is referred to again in Section |7|. 
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Figure 3. Comparison bet ween the blue profil e prese nted in this 
work with those given by 'K illeen fc Bick nell* Jl98i ' 



(open cir- 



3 rayc 

cles). ICaon et al.l ||T994) (filled circles) and Hilke r et al.l J1999 
(squares). This last authors V profile was transformed adopting 
a constant {B - V) = 0.93. 



Figure 4. Areal density of the NGC 1399 red globular clusters. 
Open circles from lDirch et all i20()j) data. Fill ed circles from the 
(scaled) density profile given by TpWhes et alj tl998). The solid 
line is the approximate fit of a Hubble law with rc = 60 arcsec 
and a central density of 72 globulars per sq. arcmin (down to 
magnitude R=23.2). 



5 THE AREAL DENSITY PROFILES 

Numerous studies of GCS in luminous elliptical galaxies 
agree in the sense that the behavior of the areal densities 
of the the blue and red GC show a similar behaviour from 
galaxy to galaxy: blue clusters usually show a more shal- 
low distribution than that corresponding to the red ones 
which, in turn, show density slopes more similar (but not al- 
ways identical) to that of the galaxy haloes. A custommary 
approach uses r^''^ laws to fit the density run with galac- 
tocentric radius although these approximations fail in the 
inner regions where GCS exhibit flat cores (see for example, 
iKundu fc Whitmordll99ilLee. Lee fc Gibsonl Eo02). These 
core profiles show a marked difference with the galaxy light 
profiles. This last feature has been used to estimate the rate 
of GC destruction due t o gravitational effects in NGC 1399 
and in o ther galaxies fCapuzzo-Do lcetta fc Donnarummal 
l200lllCap uzzo-Dolcctta fc Tcssori 199^! 

Figure|l|depicts the areal density of the red GC derived 
from D2003 data. An attempt to derive the ellipticity of the 
red clusters distribution (and also for the blue ones) based 
on azimuthal counts, indicates a value close to zero (com- 
parable with that of the galaxy itself), which we adopted in 
what follows. The inner values (filled circles) were derived 
from HST observations given bv lForbes et alJ lll998l) . Their 
counts were normalized to a limiting magnitude Ti — 23.2 
and scaled down by 2/3. This last factor takes into account 
that only this fraction of the GC seem to belong to the red 
population on the basis of the colour statistics in the in- 
ner region o f the galaxy, also base d on HST observations, 
presented bv lCrillmair 

For comparison purposes only. Figure 21 also shows a 



Hubble law characterized by a central density of 72 GC per 
sq. arcmin and a core radius of 60 arcsec. 

The blue GC density profile is depicted in Figure|^that 
also includes the best formal fitting with an r^''* law. In this 
case we note that there is a marginal trend in the sense 
that the bluest clusters((C — Ti) = 0.90 to 1.20, also shown 
in the same figure) exhibit a shallower core. As a reference 
we tentatively fit an isothermal profile iBinnev fc Tremainj 
E987) with a core radius of 275 arcsec for these objects. 
However this feature will require further confirmation using 
a deeper sampling and, in what follows, we consider the blue 
GC as a single family. 

One important aspect related with the determination 
of the GC density profiles at large galactocenric radii is the 
impact of unresolved field interlopers with colours similar 
to those of the clusters. We stress that D2003 data were 
processed with a highly successful procedure for the iden- 
tification of non-resolved sources, as the GC are expected 
to be at the Fornax cluster distance. Further spectroscopic 
data presented in R2004 indicate a contamination level on 
the order of 10 % or less. In addition, their comparison field 
(30 arcmin on a side) located three degrees towards the NE 
of NGC 1399, provides a very solid estimate of the back- 
ground level in terms of the uncertainties of the statistical 
counts (see Section 7). 
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Figure 5. Areal density of the NGC 1399 blue globular clusters. 
The straight line is the best formal fitting with an r^^^ law. The 
lower set of values corresponds to the bluest globulars (see text) 
fit with an isothermal profile with a core radius of 275 arcsec. 

6 THE GC VOLUMETRIC DENSITY 
PROFILES 

Volumetric density profiles allow an estimate of the total 
mass within a given spatial galactocentric radius as well as 
a comparison with other galaxy components (e.g. hot gas, 
dark matter) and, also, to estimate the effect of an even- 
tual (spatial) colour gradient on the projected integrated 
colours. We note that neither the galaxy light profile nor 
the GC areal densities can be approximated by means of 
single (spatial) power laws {r") except for limited ranges in 
galactocentric radius. In this work, and following R2004, we 
adopt integrable volumetric density functions of the form: 

Dens{r) = p.{riy ' (6) 

and aim to determine the and parameters that, once 
Dens(r) is integrated along the line of sight, give a good rep- 
resentation of the projected (or areal) density profiles. The 
integration was carried out within a limiting radius equal 
to the virial radius derived by R2004, Rvir ~ 5800 arcsec. 
We note, however, that changing this limit by a consider- 
able amount does not have a large effect on the results of 
the integration. 

The best formal fits to the areal density of the red GC 
are displayed in Figure |^ The results show that different 
combinations of the C—Va parameters can yield almost equiv- 
alent fits within the range of the fit (120 to 1200 arcsec) and 
that this degeneracy breaks up only in the innermost re- 
gions of the galaxy. However, deciding which of these profiles 
might be the best representation of the original GC distri- 
bution is probably impossible as the disrupting gravitational 
effects are expected to be important in these regions. 
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Figure 6. Projected volumetric density profiles that give an ad- 
equate fit to the observed red GC surface density (circles) for a 
given combination of and ^ values (see text); (J = 1.5, = 70 
arcsec (short dashed); = 1.0, = 50 arcsec (long dashed); 
C = 0.0, = 30 arcsec (solid line). The profiles overlap for galac- 
tocentric radii larger than 80 arcsec. 
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Figure 7. Same as Figure 6 but for the blue GC (open circles), 
f = 0.0, rs = 200 arcsec (continuous line); f = 1.0, = 400 
arcsec (long dashed); f = 1.5, = 700 arcsec (short dashed). 
Profiles overlap for galactocentric radii larger than 150 arcsec. 
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Figure |7| shows the same kind of diagram for the blue 



GC. 



7 RESULTS 

This section presents the results of fitting the composite 
stellar population profile described in Section |3] to the blue 
surface brightness presented in Section |1] The quality of the 
fit is estimated in terms of: 



Shape of the surface brightness profile 
Colour gradient 

Cumulative GC specific frequency 



7.1 The shape of the surface brightness profile 

In order to derive the best fitting to the observed blue sur- 
face brightness profile (defined as the one that minimizes the 
square of the brightness residuals) we proceeded as folows: 

- A grid of volumetric density profiles that, once pro- 
jected, provide acceptable fits (with residuals within 1.5 
times the mimimum dispersion) to the observed areal den- 
sities was generated for each GC family (setting =0.0, 1.0 
and 1.5) and changing the parameter. 

- These profiles were projected on the sky (from r=0 to 
1200 arcsec) and then combined within the equation (de- 
scribed in Section I^J. The constant C{B) was iteratively 
changed until the quantity k — /ib + 2.5log[N{RGC') + 
N{BGC)/C{B)\ remained constant with galactocentric ra- 
dius, yielding a value of 33.70. The mimimum dispersion of 
this quantity (±0.035) was obtained by combining GC den- 
sity profiles characterized by = 1.5 and Ts = 30 arcsec for 
the red GCS, and ^ = 1.0 and = 375 arcsec for the blue 
GC. 

- The k value was equated to [V — Mv)o + A{B) + 
2.5log[SB{RGC)] and, by adopting the corrected distance 
modulus and interstellar extinction discussed in Section 
the instrinsic specific frequency of the red GC was derived. 
The specific frequency of the blue clusters follows from C{B) 
and Sb{RGC). 

This procedure leads to C(B) = 3.5 ±0.5, Sb{RGC) = 
7.8, Sb{BGC) = 27.3. The uncertainty for the C{B) pa- 
rameter comes, mainly, from the range of values allowed by 
the errors in the outermost region of the brightness profile. 
The resulting brightness profile is compared with the obser- 
vations in Figure 8 and the fit residuals are shown in Figure 
9. This last diagram indicates that, from 30 to 720 arcsec 
the approach gives a very good fit to the observations with 
maximum deviations of about ±0.05 mag and an overal rms 
of ±0.035 mag. 

The inner 30 arcsec show a deviation in the sense that 
the nuclear region is brighter than the model. Improving the 
fit in this region might require a higher exponent (and also 
the consideration of seeing effects). However, as discussed 
later, some arguments support the presence of a distinct 
stellar population in the galaxy nucleus. 

The volumetric profiles assumed a spherical geometry. 
Even though the nature of the residuals in Figure 9 cannot 
be completely asserted, we note that they seem correlated 
with the ellipticity variations depicted in the upper curve 




400 

(arcsec) 



Figure 8. Best profile fit to the observed blue surface brightness 
of NGC 1399 (solid line) The short dashed line corresponds to the 
light contribution of the "red" stellar population; the long dashed 
line is the light contribution from the "blue" stellar population. 
The associated GC have S„{RGC) = 3.3 and S„iBGC) = 14.3 



of that figure. The assumption of an oblate geometry did 
not improve the fit suggesting the necessity of a three axis 
ellipsoid for a further comparison. 

We stress that, due to the small fractional contribution 
of the blue stellar population to the integrated light in the 
inner regions of the galaxy, the fit is rather insensitive to the 
set of rs — C values adopted for the blue GC. The reason for 
adopting the quoted values are commented on in Section 11. 



7.2 Sources of uncertainty 

One source of uncertainty comes from the background level 
adopted as representative of the field contamination by non- 
resolved objects with colours in the domain of the GC 
colours. The comparison field used by D2003 indicates val- 
ues of 0.31 and 0.08 objects per sq. arcmin for the blue and 
red GC respectively. This field (approx. 35 by 35 arcmin) 
includes a large number of objects and yields low counting 
uncertainties (±0.016 and ±0.008). However, the doubt re- 
mains in the sense if this is in fact representative of that of 
the central regions of the Fornax cluster (three degrees away 
to the SW). 

An increase (decrease) of the background levels by 50% 
of the adopted values changes the shape of the density (and 
brightness) profiles forcing an increase (decrease) of the 
C{B) ratio of ~ 0.5 in order to keep an acceptable fit to 
the observed profile. 

The combined uncertainty of the quantities involved 
(fc : ±0.035); A{B) : ±0.02; (K-A'/„) ; ±0.1) lead to ±0.8 in 
the Sb (RGG) parameter. This value, and the quoted uncer- 
tainty on the G{B) ratio, in turn yield ±4.8 for Sb{BGG). 
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Figure 9. Blue profile fit residuals as a function of the elliptical 
galactocentric radius (solid line). The maximum deviations 
are within ±0.05 mag with an overall rms of ±0.035 mag over 
a galactocentric radius from 30 to 730 arcsec. Note the marked 
deviation within 30 arcsec. The upper (dashed) curve shows the 
variations of the galaxy fiattening (q = b/a). A constant q value 
(0.83) was adopted for re larger than 420 arcsec. 



7.3 [B — Rkc) colour gradient 

A composite red profile was computed by adding tlie red 
liglit contribution of eacli stellar population (adopting the 
{B — Rkc) colour indices and colour excess discussed in 
Section Tlie resulting {B — Rkc) colour gradient is dis- 
played in Figure ITHl that also includes the observed values 
listed in Table 1. This diagram shows that, at galactocen- 
tric radii larger than 30 arcsec, the predicted colours are 
consistent with the observations. In turn, these ob servations 
are in very good agreement with iMichardI i200[ l). This last 
author gives, at a galactocentric radius ro = 22.8 arcsec, 
{B — Rkc) ~ 1-54 (0.02 mag redder than ours) and a colour 
gradient dlog(B — Rkc) / dlogij) — —0.05 that compares well 
both with the observed (— 0.06±0.02) and with the predicted 
(-0.045) colour gradient. 

Another estimate of the {B — Rkc) gradient was ob- 
tained by combining our blue profile and the Rkc one given 
in D2003. These last authors adopted a "preferred" sky level 
and warned about the high sensitivity of the profile at the 
outermost radius. The uncertainty of their profile becomes 
~ ±0.05 mag within a galactocentric radius of 8 arcmin that 
we adopt as the boundary of our combined {B — Rkc) colour 
profile. The resulting {B — Rkc) is also depicted in Figure lTOl 
after applying a zero point shift of -0.1 mag. This correc- 
tion is necessar y in order to bri ng the resulting colours into 
agreement with lMichardI ll2000l) and with ours. 

As a summary, both our observed [B — Rkc) gradient 
and that derived by combining our blue profile with the red 




(arcsec) 



Figure 10. Galaxy [B — Rkc) colour as a function of elliptical 
galactocentric radius Tc (solid curve). The straight line is the ex- 
pected gradient from the model described in S ection 2. The squ are 
indicates the {B — Rkc) colour measured bv iMichardI i2OO0t) at 
ro = 22.8 arcsec. The bar at 228.8 arcsec indicates the uncer- 
tainty of our colour colour gradient. Open circles show the result- 
ing colour gradient obtained by combining our photometry and 
that given in D2003 (see text). 



one from D2003, are consistent with that predicted by the 
profile fit. 



7.4 The intrinsic and cumulative GC specific 
frequencies 

The intrinsic specific frequencies derived for each GC family 
were transformed to the V band, adopting the [B — V)o 
colours given in Section El yielding Sn{RGC) = 3.3 ± 
0.3 and Sn{BGC) = 14.3 ± 2.5. In the case of the 
red GC, the derived specific frequency is coincident with 
the representative value quoted for normal ellipticals (see 
iHarris. Harris &i McLaughlin Ill99d) . The composite cumu- 
lative (with galactocentric radius) specific frequency was 
computed adopting a total number of 700 GC within a 
galactoc entric radius of 120 arcsec (from the density curve 
given bv lForbes et al.lll99d) . GC counts at larger radii were 
derived fom D2003 data adopting a limiting magnitude 
Ti — 23.2, a fully gaussian integrated luminosity function 
with a turn-over at this magnitude, and an overall complete- 
ness of 90 %. In turn, the galaxy V surface brightness profile 
was derived from our observed blue profile and the {B — V) 
colour gradient computed from the model fit. The adoption 
of the distance modulus and interstellar extinction discussed 
in Section 121 lead to the Sn values shown in Figure [TTI The 
derived cumulative Sn is somewhat smaller, but comparable 
within the uncer tainties, with results given in D2003 a nd in 
previous works JOstrov et alJll998l : iForte et al.ll2002h that 
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Figure 11. Cumulative specific frequency of the GCS associated 
witli NGC 1399. Open circles represent Sn values derived from 
the actual counts (transformed to total GC population; see text) 
combined with the V surface brightness profile predicted by the 
model. The solid line comes from both the fitted GC areal densi- 
ties and the model V profile. 

have pointed out that the NGC 1399 GCS cannot be con- 
sidered as an anomalously high 5„ system. 



8 A CAVEAT ABOUT THE APPROACH 

An alternative fitting to the observed surface brightness pro- 
file can be obtained by assuming that the stellar halo follows 
a density profile similar to that of the red GC and forcing a 
spatial colour gradient d{B — Rkc) / d{logr) = —0.045. 

This blueward trend with increasing galactocentric ra- 
dius, if due to a change of metallicity, would increase the B 
brightness of the stellar population (at constant mass) in the 
outer reg i ons. T he luminosity increase was estimated from 
IWorthevI il994l) population synthesis models and assuming 
an age of 12 to 17 Gy (the results are not very sensitive 
to age in this regime) . The forced spatial gradient also pro- 
vides a good fit to the observed projected {B — Rkc) colour 
gradient. 

If that gradient were also shared by the red GC, their 
(projected on the sky) number vs. colour statistic, computed 
on the basis of the volumetric density profile, would show 
a modal colour shift of -0.14 mag in (B — Rkc), and close 
to -0.20 mag in (C — Ti), over a galactocentric radius of 20 
arcmin. This result is in confiict with the observations (e.g. 
D2003) that show a practically constant modal colour with 
galactocentric radius for both GC populations. 

Despite this last comment, the main caveat about the 
profile fit approach presented in this work is the possibility 
that the chemical enrichment history of the red field stellar 
population might be in fact decoupled from that of the red 



GC. In fact there are theoretical reasons to expect some de- 
gree of decoupling. For example, if GC form at the beguining 
of the starburst when external pressure and content of gas is 
highest; as the star formation rate (SFR) decreases, chem- 
ical enrichment may continue to ocur (but no new GC are 
formed). This would presumably lead to a color difference 
between GC and the galaxy halo. We note, however, that in 
the case of the halo, the result will be luminosity weighted, 
i.e., the effect on the color will depend upon the ratio of 
the number of stars formed during the SFR peak and those 
formed subsequently. If this ratio is high, the resulting halo 
color could be only slightly redder than those of the RGC. 



9 THE STELLAR METALLICITY 
DISTRIBUTION FUNCTION 

This section presents a tentative picture about the expected 
stellar metallicity distribution function (MDF) based on the 
two component model presented above. Observations of re- 
solved stars in the halo of NGC 5128 show marked differ- 
ences of their metallicity distribution whe n compared with 
that of GC in number by number statistics iHarris fc HarrisI 
l2002h . 

If in fact the GC specific frequency is the result of a 
composite population, as suggested by the previous section, 
an estimate of the number ratio of metal rich (red popula- 
tion) to metal poor (blue population) stars can be derived 
through the intrinsic specific frequencies and their mass to 
light ratios. 

The expected stellar mass associated with A'^ clus- 
ters within a (projected or spatial) galactocentric radius is 
M{stellar) = {M/L){N/S„). For the red p opulation we 
adopt a blue {M/L)redpop = 10 derived bv ISaglia et alJ 
|2000) for the inner regions of NGC 1399 (where that pop- 
ulation is dominant) on the basis of their kinematics study. 
By assuming that the colour difference between both stellar 
populations arises, mostly, from a metallicity difference close 
to 1.0 dex (see below)j and the results from stellar popula- 
tion modelling (e.g. IWorthev.1994'1 . we get {M/L)biue pop = 
5.75 for the blue population. This ratio is consitent with the 
fact that the turn over of the integrated luminosity function 
of the blue GC in the blue band is about .4 mag brighter 
than t hat of the red GC as determined bv lGrillmair et all 

The resulting stellar mass ratio then is: 

M{Red) _ Sb{BGC) {M/L)rcdpop NjRGC) 

M{Blue) Sb{RGC){M/L)uu.popN{BGC) ^' 

In the particular case of the region between 120 and 420 arc- 
sec of galactocentric radius (see Figure^ where N{RGC) ~ 
N{BGC), the adoption of the specific frequencies derived 
in Section |7| leads to a mass ratio close to 6. If the stel- 
lar mass distribution function of both populations do not 
differ markedly, this result implies that blue stars will be 
relatively inconspicuous in a number statistic that will be 
largely dominated by the red population stars. 

The shape of the MDF for these stars can the be es- 
timated by transforming the GC colours to metallicity and 
combining both populations according to the mass scaling 
factor derived above. This procedure deserves a cautious 
comment taking into account several factors that may have 
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Figure 12. Expected metallicity distribution function (MDF) 
for the stellar population within a (projected) galactocentric ra- 
dius of 2 to 7 arcmin whose colour distribution is shown in Fig- 
ure O The predicted composite MDF shows a good qualitative 
agreement with observ ations of resolved halo stars in NGC 5128 
iHarris fc Harri£l2002l) . i.e. a dominant and wide high metallicity 
population and a less conspicuous and extended blue tail. 



an impact on colour-metallicity relat ions, like the morp hol- 
ogy of the horizontal branch stars iGeisler et aljll996ll or 
the existence of an age spectrum among the NGC 1399 GC 
(Forbes et al. 2001). 

Integrated (C — Ti) colours for the GC sample shown 
in Figure that covers a galactocentric radius from 11 to 
38 Kpc, were transformed to metallicity by means of the 
calibration presented bv .Harris fc Harris; (■2002i'l 

[Fe/H] = -6.037[l.0 - 0.82(C - Ti)o + 0.162(C - (8) 

that improves the linear relation given by ICeisler fc Fortd 
il99Qil . The GC sample was divided in two groups: blue GC 
with (C-Ti) = 0.9 to 1.60 and red GC with (C-Ti) = 1.5 
to 2.2 (i.e. allowing a 0.1 mag overlap between both samples) 
and correcting for interstellar reddening. 

The resulting MDF (in relative units), after applying 
the scaling factor described before, is shown in Figure 1121 
This distribution shows a good qualitative agreement with 
that corresponding to the outer fields in NGC 5128 from 
[Har ris fc Harris (2002), i.e, a dominant red population with 
a tail extending to low metallicity values. The results of these 
last authors stand on color magnitude diagrams of resolved 
stars interpreted in terms of evolutionary tracks with differ- 
ent ages and metallicities. In turn, the MDF presented in 
this work relies on an empirical integrated-colors vs. metal- 
licity relation. As a consequence, any further quantitative 
comparison will require a previous effort in order to check 
for the consistency of both approachs in terms, for example, 
of the metallicity scale. 



10 STELLAR MASS AND GC FORMATION 
EFFICIENCY 

The GC formation efficiency is a key parameter in the un- 
derstanding of the link between these clusters an the diffuse 
stellar population of the halo. For a given GC specific fre- 
quency, the formation efficiency, in terms of the associated 
stellar mass, is 



{M/L)c 



(9) 



where (M/L) is the mass to luminosity ratio of the pop- 
ulation, Mgc is the mean GC mass (2.4 X lO'' Mq), and 
the numerical factor c is the luminosity corresponding to 
M\ = —15.0 in solar luminosity units. 

Adopting the specific frequencies derived in Section 
7, the mass to luminosity ratios discussed in the previ- 
ous section, with c = 1.55 x 10* for the blue band, yield 
mac = 1.25 X 10"^ and tibgc = 7.60 x 10" ^ 

Beside, and following lMcLaughlin|[T999l this definition 
can be revised in terms of the total baryonic mass, i.e., in- 
cluding the mass of hot gas. In this case we define a shell 
with inner and outer (spatial) radii of 11 and 66 Kpc re- 
spectively. The adoption of the inner boundary attemps to 
decrease the effect of gravitational disruption on the GC 
sample while the outer one is the limit of our brightness 
profile. About 1800 RGC and 1500 BGC fall within these 
boundaries that enclose 3.54 x lO" and 0.55 x 10^^ Mq of 
stars of the red and blue populations, respectively. On the 
other side, adopting the X ray analysis by Jones et al. (1997) 
and correcting their adopted distance (24 Mpc) to ours, we 
obtain a hot gas mass of 0.28 x IO^^Mq. These quantities 
lead to an overall efficiency r) — 1.8 x 10~^ for both GC 
populations combined. This value is somewhat smaller, but 
compares well, with the GC formation efficiency (2.9 x 10""^) 
derived by .McLaughlin, (.1999.1 . 
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A COMPARISON BETWEEN GLOBULAR 
CLUSTERS, STARS, HOT GAS AND DARK 
MATTER DISTRIBUTIONS 



11.1 Red Globular Clusters and Total Galaxy 
Mass 

The results presented in the previous sections are consis- 
tent with the existence of a link between the RGC and the 
(red and dominant) halo stellar population of NGC 1399. 
This connection is also supported by the GC kinematic be- 
haviour. For example, ISaelia et alJ i2000fl give a radial ve- 
locity dispersion for the halo stars of about 280 km 
(at 120 arcsec from the centre) in good agreement with the 
RGC value of 265 km (from R2004 data, without the 
velocity cut-off adopted by these authors). Taking into ac- 
count that the derived volumetric profile for these clusters 
is consistent with a logarithmic spatial density derivative 
{dln{p) / dlnr) = —2.9 (at galactocentric radii larger than 
~ 100 arcsec), and using the Jeans equation for the spheri- 
cal and isothermal case, the resulting total mass within 100 
Kpc of galactocentric radius ranges from 4.4 to 4.8 x W^^ Mq 
(for an anisotropy parameter (3 from 1.0 to —0.5). 
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Figure 13. Upper panel: The NGC 1399 blue globular clusters 
density profile (open circles) compare d with that of the hot X 
ray emitt ing gas (continuou s line: from lPaolillo et~ .120021: dahed 
line: from ljones et alll997l . Lower panel: blue GC compared with 
the (projected) dark matter profile obtained by R2004 using GC 
kinematics. 



sets oi ( — rs values, we adopted ( = 1.0 and rs = 375 arcsec 
as adequate for the BGC noting that our scale length can 
be considered as identical (to within the uncertainties) with 
the value given in R2004. 

As discussed in Section |n] volumetric profiles with dif- 
ferent values of the parameter can yield almost equivalent 
fits to the observed areal densities for galactocentric radii 
larger than 2 arcmin. This means that, beyond the very good 
agreement between the distribution of the blue GC and dark 
matter, the data do not allow us to support unambiguously 
that an NFW profile also holds in the inner regions of the 
galaxy. However, all these profiles are consistent with an 
spatial derivative dln{p) / dln{r) ~ —1.4 from 100 arcsec and 
outwards in galactocentric radius. If this value is adopted as 
representaive for the hot gas (according to Figure 13) and 
for a temperature range of 1.1 to 1.4 Kev (Jones at al. 1997; 
Paolillo et al. 2002), the total gravitating mass within 100 
Kpc assumig hydrostatic gas equilibrium, becomes 5.3 to 
5.8 X IO^^Mq, i.e., comparable with that derived from the 
RGC kinematics mentioned before. 

The fact that the slope of the density profiles of the 
blue GC and that corresponding to the X ray emission asso- 
cia ted with NGC 139 9 seem similar was already suggested 
bv lHilker et alJ 1^23) ■ turn, the similarity between the 
distribution of the blue GC and that of dark matter was 
also pointed out by ICote et alj il998l) in their analysis of 
the outer regions of the M87 GCS. 



11.2 Blue Globular Clusters, Hot Gas and Dark 
Matter 

The X ray emission of t he Fornax cluster has been the sub- 
ject o f several works (e.g. lSerlemitsos et al.ll99Sllkebe et alJ 
[1992). ROSAT observations have been discussed by Jones et 
al. (1997) who fit a surface brightness profile of the form: 



S{r) = 5(0) 



1 + 



-1 -3/3+0.5 



(10) 



centered on NGC 1399, obtaining be st fitting parameters 
a = 45 arcsec and /3 = 0.35. In turn, IPaolillo et aP ||2002|) 
give a thorough description of the complex X ray emission 
and suggest a composite beta model (includin g "central" , 
"galax y" a nd "cluster" components). Both the ljones et alJ 
^1997^ and lPaolillo et a l. ( 2002) profiles are depicted in Fig- 
ure 1131 (upper panel) and exhibit a good coincidence with 
the surface density of the blue GC. 

The lower panel in Figure 13 shows the projected pro- 
file that gives a representation of the dark matter dis- 
tribution (also shifted in ordinates) derived by R2004 on 
the basis of their GC kinematic study . These last authors 
note that a NFW profile (i.e. C = 1.0: iNavarro et aljlT99^ 
iNavarro. Frenk fc White Il997^ with a scale length = 394 
arccsec gives a good representation of the inferred volumet- 
ric density of dark matter. We remark that, although their 
fit covers a galactocentric range from approximately 2 to 8 
arcmin (re^^"* Irom 3.31 to 4.68), the profile extrapolation 
holds up to 20 arcmin (i.e. the approximate limit of the 
D2003 data). For this reason, and among all the acceptable 



12 DISCUSSION 

The surface brightness profiles presented in Section |1] can- 
not be matched with a single r^^'^ law, a feat ure already 
noted in previous works (e.g. ISchombertlll986i) that have 
also pointed out the presence of an outer corona that would 
entitle the classification of NGC 1399 as a cD galaxy (but 
see the difficulties inherent to such a classification described 
in R2004). 

If in fact the formation of a given GC family occurs 
concomitantly with that of a diffuse stellar population that 
shares its age and chemical abundance, it may be expected 
that, after gravitational relaxation, both stars and clusters 
will exhibit similar spatial distributions. Although dynami- 
cal friction might play a differential role that remains to be 
estimated, globulars could be identified as tracers of both 
the spatial distribution and of the kinematic behaviour of 
coeval diffuse stellar populations. Then, the galaxy bright- 
ness profile would be the result of the "luminosity weighted" 
composition of those events. This approach gives a proper 
fit of the NGC 1399 profile shape by assuming the presence 
of two dominant populations (although they might exhibit 
their own degree of internal inhomogeneity in terms of age 
and chemical abundance) typified by very different intrinsic 
specific frequencies 5'„. As the blue stellar population be- 
comes more prominent outwards, due to its shallower areal 
density slope, the degree of "cDness" of the galaxy can also 
be identifyied as a result of that population. 

An important result of the analysis is that, even though 
both the red and blue GC families do not exhibit a de- 
tectable change of their modal colours over a galactocentric 
range close to 100 Kpc, the composition of their associated 
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diffuse stellar populations do lead to a blueward colour gra- 
dient that compares well with the observations. 

A caveat remains, however, in the sense that the galaxy 
profile can be matched by adopting a stellar population den- 
sity profile similar to that of the red GC and forcing a spatial 
colour gradient. This would imply, on the one hand, that for 
some reason the stellar and red GC enrichments were decou- 
pled and, on the other, that blue GC formed in a process 
unable to form non-clustered stars. 

Although the overall globular cluster formation effi- 
ciency obtained in previous sections seems comparable to 
other values available in the literature (e.g. McLaughl^ 
^9^) the considerably higher intrinsic S„ value of the blue 
population in comparison with the red one also implies a 
higher globular formation efficiency that may be indicative 
of a distinct formation environment. Possible explanations 
of this situation have been given in the literature in con- 
nection with the role of re-ionization as responsible of the 
ter mination of the star f orming process or in the frame of 
the ISearle fc ZinnI lll978f) scenario. In this last case, rela- 
tively small gas clumps associated with the formation of the 
BGC are disrupted by Supernovae. 

As noted before, the profile fit does not give an ac- 
ceptable description of the inner 30 arcsec in galactocentric 
radius where the galaxy is brighter and redder than pre- 
dicted. This situation may arise from the use of an inappro- 
priate analytical profile, and/or the effect o f GC disruption 
by g ravitational forces (see, for example, IVesperini et alJ 
l 2003r. Apert u re int egrated (C — Ti) colours measured by 
lOstrov et aP are indicati ve of the presence of a 

high mctallicity stellar population JPickle!jl985l : lKuntschneil 
[2OOQ ). This si tuation has similarities w ith other galaxies 
(e.g. the MW; iMcWiUiam fc RichllTooll where bulge stars 
seem to have metallicities somewhat higher than those of 
the RGC. A distinctive kinematic behavior of the NGC 1399 
central regio n can be seen in th e velocity dispersion profile 
obtained bv ISaelia et alJ (|200(J) . These facts indicate that 
the galaxy nucleus may have harbored a different sequence 
of stellar formation events disparate from those observed 
over much larger spatial scale lengths. 



13 CONCLUSIONS: A TENTATIVE SCENARIO 

Previous sections have described a possible approach to the 
link between GC and the diffuse stellar population in NGC 
1399 on a purely photomet ric basis. Although spectroscopic 
results jForbes et al.ll200 J) indicate the existence of an age 
spectrum among its GC (and, presumably, in the associated 
diffuse stellar populations), the analysis of a significant vol- 
ume of spectroscopic data in terms of age and metallicity is 
still missing. 

In a cautious way, however, the results presented in the 
previous sections, and other available in the literature, can 
be ordered into a working scenario that attempts to delin- 
eate a temporal sequence of events that might end in the 
observed situation, namely: 

- The blue GC in NGC 1399 exhibit modal integrated 
colours that are comparable to those observed in other 
galaxies and are consistent with a frequently mentioned fact: 
The relative homogeneity of the blue GC population from 
one system to another, suggesting that their origin might 



have occurred before the formation of the dominant stel- 
l ar population of the galaxy ha lo es. However, recent works 
(Lotz. Mil ler fc Ferguson! 120041: IStrader. Brodie fc Forbes! 
j2004) indicate that a correlation between the integrated 
BGC colors and the galaxy absolute magnitude may exist. 
These results imply that, even though the BGC might have 
formed in a very early stage, these clusters "knew" about 
the galaxy they were forming in. 

- There is a remarkable coincidence of the spatial distribu- 
tion of the blue GC with that inferred for dark matter by 
R2004 on the basis of their kinematic study. We stress that, 
even though these last authors perform their analysis in a 
galactocentric range from 2 to 8 arcmin, the extrapolation 
of their dark mass profile (projected on the sky) matches 
that of the BGCs up to 20 arcmin (the outer limit in D2003 
and in our analysis). 

- X ray emitting hot gas, usually assumed to be at the virial 
temperature of galaxy clusters, also show a spatial distribu- 
tion that, within the uncertainties, is coincindent with that 
of the blue GC. We note that the m ean metallicity of this gas 
{[Fe/H] = -0.4: IJones e'ni]ll997ri is comparable to that of 
the RGC. 

- A revision of the R2QQ4 kinematic data suggests that at 
least the outer blue GC in NGC 1399 may exhibit, as in other 
galaxies, a detectable amount of "solid-body" -like rotation. 
The presence of angular momen tum in these clusters has 
been noted, for examp le, in M87 Kissler-Patig fc Gebhard3 
Il998h and NGC 4472 JCote et alJl2003!) . 

All these features strongly suggest that blue GC may 
have had their origin at a very early phase of stellar for- 
mation, sharing the evolution of dark matter clumps in the 
non dissipative collapse that finally formed the most mas- 
sive component of galaxy haloes, a fact that might also be 
reflected by their high radial velocity dispersion. In this con- 
text, the usefulness of simple closed box-typ e chemical evolu- 
tion models has been noted, for example, bv lHarris fc Harris 
j2000!) . lDurrell. Harris fc Pritchct' (2001*) and, in the partic- 
ular case of the MW, by[VanDalfson & Harris ( 2004). Along 
these lines, the scenario described by Harris fc Harris (2002) 
for NGC 5128, where early star formation goes on with an 
initial stage of rapid infall of metal poor gas, and based on 
the similarities with NGC 1399 noted in Section 9, may also 
seem relevant for this last galaxy. 

In turn, red GC seem connected with the formation of 
the dominant stellar population, probably at the expense of 
gas that has suffered a rather inhomogeneous dissipative col- 
lapse leading to angular momentum loss and to their broad 
chemical abundance distribution. This process might have 
been boosted by the effect of the very early enrichment pro- 
vided by the blue GC (and the associated diffuse stellar 
population) . 

The idea that GC may have played a role as early 
"seeds" of chemical enrichment w as already put forward by 
iDe Young. Lind fc StromI il983h . The main problem with 
their model is that the small mass locked in GC could not 
provide enough heavy elements to increase the overal metal- 
licity level. This problem is overcome if, in fact, most of the 
enrichment comes not from the cluster themselves but from a 
much more massive (although less conspicuous) diffuse stel- 
lar population. The role of the oldest GC and its associated 
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stellar population in the re-ionization of the Universe would 
also be enhanced on the same argument. 

As a reference, we note that such a blue stellar pop- 
ulation might provide some 1.9 x IO^Mq of heavy ele- 
ments within a galactocentric radius of 100 Kpc. This mass 
of heavy elements arises in an estimated stellar mass of 
8.1 X W^^Mq contained in the blue population (from the 
blue stellar light profile and the (M/L) ratios given in Sec- 
tio n 9) and adopting the net yield s (r; = 0.023) given 
by iPortinari. Chiosi fc Bressaiil lll998fl and revised by Cora 
(2004 in prep.) for a very low metallicity population. If the 
current red stellar population (8.1 x 10^'' M0) was mostly in 
gaseous form at that stage, the injected heavy metals would 
have rised the chemical abundance to about 0.15^©. 

This enrichment is about one third of the mean abun- 
dance of the red stellar population (inferred from the asso- 
ciated RGC), i.e, not enough to explain the observed dif- 
ferences between both populations. However, that level of 
enrichment might have had an important impact on the gas 
cooling function, improving the star formation efficiency and 
triggering the process that ended up with the broad abun- 
dance distribution of the red stellar population (and also on 
the current chemical composition of the hot gas). 

The much higher formation efficiency of the blue clus- 
ters suggest very different environmental conditions from 
those that prevailed for the red globulars. This fact, and the 
impact of a diffuse and low metallicity stellar population on 
the subsequent process of galaxy formation remains worthy 
of exploration. 

A final comment about the possible existence of such 
a stellar population around NGC 1399 deals with the fact 
that, being that galaxy in the neighborhood of the Fornax 
cluster bary center, it is not clear if the association is with 
this galaxy or with the cluster as a whole. For example, 
jBassino et alj|2003f) argue that GC, apparently not associ- 
ated with galaxies in the field, can be identified at angular 
distances as large as 4800 arcsecs from NGC 1399 (a distance 
comparable to the galaxy virial radius); the eventually as- 
sociated diffuse population at that radius would be hard to 
detect as the expected blue surface brightness would be in 
the range from 31 or 32 mag per sq. arcsec. 



ACKNOWLEDGMENTS 

The authors are grateful to Dr. Juergen Scheer for allow- 
ing the use of his data. One of us (FF) carried out this work 
with a Research Fellowship from CONICET . This work was 
funded with grants from CONICET, ANPCYT and Univer- 
sidad Nac. de La Plata (Argentina). 



REFERENCES 

Ashman K.M., Conti A., Zepf S.E., 1995, AJ, 110, 1164 
Bassino L., Cellone S.A., Forte J.C., Dirsch B., 2003, A&A, 
399, 489 

Baugh C. M., Cole S., Frenk C. S., Lacey C. G., 1998, ApJ, 
498, 504 

Beasley M.A., Baugh CM., Forbes D., Sharpies R.M., 
Frenk C.S.,2002, MNRAS, 347, 1150 



Binney J., Tremaine S., 1987, in Galactic Dynamics, 

Princeton Univ. Press 
Caon N., Capaccioli M., Donofrio M., 1994, A&AS, 106, 

199 

Canterna R. 1976, AJ, 81, 228 

Capuzzo-Dolcetta R., Tesseri A., 1999, MNRAS, 308, 961 
Capuzzo-Dolcetta R., Donnarumma, I., 2001, MNRAS, 
328, 645 

Cohen J. G., Blakeslee J.P., Rhyzov A., 1998, ApJ, 496, 
808 

Cote P., Marzke R.O., West M. J., 1998, ApJ, 501, 554 
Cote P., McLaughlin D.E., Cohen J.G., Blakeslee J. P. 2003, 
ApJ, 591, 850 

De Young D. S.,Lind K., Strom S.E. 1983, PASP, 95, 401 
Dirsch B., Richtler T., Geisler D., Bassino L.P., Gieren 

W.P., 2003, AJ, 125, 1908 
DurrelU P. R., Harris E. W., Pritchet C. J., 2001, AJ, 121, 

2557 

Eggen O., Lynden BeU D., Sandage A., 1962, ApJ, 136, 748 
Forbes D.A., Brodie J. P., Huchra J., 1997, AJ, 113, 887 
Forbes D.A., Brodie J.P., GriUmair C.J., 1998, MNRAS, 
293, 325 

Forbes D.A., Beasley M.A., Brodie J. P., Kissler-Patig M., 

2001, ApJ, 563, 143 
Forbes D.A., Forte J.C, 2001, MNRAS, 322, 257 
Forte J.C, Strom S.E., Strom K., 1981, ApJ, 245, L9 
Forte J.C, Geisler D., Kim E., Lee M.G., Ostrov P.G., 2002 

in Geisler D., Grebel E.K. & Minniti D.,eds., Proceedings 

of the 207th lAU Symposium, Extragalactic Star Clusters, 

Astron. Soc. Pac, San Francisco, p. 251 
Geisler D., 1996, AJ, 111, 480 
Geisler D., Forte J.C, 1990, ApJ, 350, L5 
Geisler D., Lee, M. C, Kim E. , 1996, AJ, 111, 1529 
GriUmair C.J., Forbes D.A., Brodie J.P., Elson R.A.W., 

1999, AJ, 117, 167 
Harris H., Canterna R., 1977, AJ, 82, 798 
Harris W.E., 1996, AJ, 112, 1487 
Harris W.E., van den Bergh S., 1981, AJ, 86, 1627 
Harris W.E., Harris G.L.H, McLaughlin D.E., 1998, AJ, 

115, 1801 

Harris G. L. H., Harris W. E., 2000, AJ, 120, 2423 
Harris W.E., Harris G.L.H. , 2002, A J, 123, 3108 
Harris, W. E. 2003, Extragalactic Globular Cluster Sys- 
tems, Proceedings of the ESO Workshop held in Garching, 
Germany, Agust 27-30 2002, p. 317. 
Hilker M., Infante L., Richtler T., 1999, A&AS, 138,55 
Ikebe Y. et al. 1992, ApJ, 384, L5 

Jones C, Stern C, Forman W., Breen J., David L., Tucker 

W., 1997, ApJ, 482, 143 
KiUeen N.E.B., BickneU G.V., 1988, ApJ, 325, 165 
Kissler-Patig M., Gebhardt K. 1998, AJ, 116, 2237 
Kissler-Patig, M. 2002, in Geisler, D., Grebel E.K, & Min- 
niti D., eds.. Proceedings of the lAU 207th Symp., Extra- 
galactic Star Clusters, Astron. Soc. Pac, San Francisco, 
p. 207 

Kissler-Patig, M. 2003, Extragalactic Globular Cluster Sys- 
tems, Proceedings of the ESO Workshop held in Garching, 
Germany, Agust 27-30 2002, ESA. 
Kundu A., Whitmore, B.C., 1998, AJ, 116, 2841 
Kuntschner H., 2000, MNRAS, 315, 184 
Larsen S. S., Richtler T., 2000, A&A, 354, 836L 
Lee H-c, Lee Y.W., Gibson B.K., 2002, AJ, 124, 2664 



Globular Clusters as Tracers of Stellar Bi-modality in Elliptical Galaxies 



Lotz J. M., Miller B. W., Ferguson H. C, 2004, AJ, 613, 
262 

McLaughlin D.E., 1999, AJ, 117, 2398 
McWilliam, A., Rich, R. M., 1994, ApJS, 91, 749 
Michard R., 2000, A&A, 360, 85 

Navarro J.F., Frank C.S., White S.D.M., 1996, ApJ, 462, 
563 

Navarro J.F., Frank C.S., White S.D.M., 1997, ApJ, 490, 
493 

Ostrov RG., Gcislcr D., Forte J.C., 1993, AJ, 105, 1762 
Ostrov P.G., Forte J.C., Geisler D., 1998, AJ, 116, 2854 
Paolillo M., Fabbiano G., Peres G., Kim D.-W. 2002, ApJ, 
565, 883 

Peebles P.J.E., Dickc R.H. 1968, ApJ, 154, 891 
Pickles A. J., 1985, ApJ, 296, 340 

Portinari L., Chiosi C., Bressan A., 1998, A&A, 334, 505 
Reed B.C., Hesser J.E., Shawl S.J., 1988, PASP, 100, 545 
Richtlcr T., Dirsch B., Geislor D., Gcbhardt K. Hilker M., 
Alonso M.V., Forte J.C., Grobel E.K., Infante L., Larsen 
S., Minniti D., Rojkuba M., 2004, AJ, 127, 2094 
Ricotti, M. 2002, MNRAS, 336, L33 

Saglia R.P., Kronawitter, A., Gerhard, O., Bender, R., 
2000, AJ, 119, 153 

Santos M.R. 2003, in Extragalactic Globular Cluster Sys- 
tems, Proceedings of the ESO Workshop held in Germany, 
27-30 Aug. 2002, p. 348 

Schlegel D., Finkbeiner D., Davis M., 1998, ApJ, 500, 525 

Schccr J. 1987, Appl. Optics, 26, 3077 

Schombert J.M., 1986, ApJS, 60, 603 

Schweizer F., 1987, Nearby Normal Galaxies. From the 
Planck Time to the Present. Springer, New York, p. 18. 

Searle L., Zinn R., 1978, ApJ, 205, 357 

Sorlcrnitsos P., Locwenstcin M., Mushotzky R., Marshall 
F., Petre R., 1993, ApJ, 413, 518 

Strader J., Brodie J. P., Forbes D. A., 2004, AJ, 127, 3431 

VanDalfscn M. L., Harris E. W., 2004, AJ, 127, 368 

Vespcrini E., Zcpf S.E., Kuudu A., Ashman K.M., 2003, 
ApJ, 593, 760 

Worthey G., 1994, ApJS, 95, 107 



